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Abstract 

Analytical solutions of the Schrodinger equation for the one-dimensional 
quantum well with all possible permutations of the Dirichlet and Neu¬ 
mann boundary conditions (BCs) in perpendicular to the interfaces 
uniform electric held S are used for the comparative investigation of 
their interaction and its influence on the properties of the system. 
Limiting cases of the weak and strong voltages allow an easy math¬ 
ematical treatment and its clear physical explanation; in particular, 
for the small S^ the perturbation theory derives for all geometries 
a linear dependence of the polarization on the held with the BC- 
dependent proportionality coefficient being positive (negative) for the 
ground (excited) states. Simple two-level approximation elementary 
explains the negative polarizations as a result of the held-induced de¬ 
structive interference of the unperturbed modes and shows that in 
this case the admixture of only the neighboring states plays a domi¬ 
nant role. Different magnitudes of the polarization for different BCs 
in this regime are explained physically and conhrmed numerically. 
Hellmann-Feynman theorem reveals a fundamental relation between 
the polarization and the speed of the energy change with the held. It 

*King Abdullah Institute for Nanotechnology, King Saud University, P.O. Box 2455, 
Riyadh 11451 Saudi Arabia; E-mail; oolendski@ksu.edu.sa 


1 



is proved that zero-voltage position entropies Sx are BC independent 
and for all states but the ground Neumann level (which has Sx = 0) 
are equal to In 2 — 1 while the momentum entropies depend on the 
edge requirements and the level. Varying electric field changes po¬ 
sition and momentum entropies in the opposite directions such that 
the entropic uncertainty relation is satisfied. Other physical quantities 
such as the BC-dependent zero-energy and zero-polarization fields are 
also studied both numerically and analytically. Applications to differ¬ 
ent branches of physics, such as ocean fluid dynamics and atmospheric 
and metallic waveguide electrodynamics, are discussed. 


1 Introduction 


One-dimensional (ID) qnantnm well (QW) with vanishing at its edges wave- 
fnnction \I^(a;) is an example that is included in every textbook on quantum 
mechanics [1]. Influence on its properties of an electric field that is applied 
perpendicularly to the confining planes [2-7] is of not only purely academic in¬ 
terest but also of the large technological applications. By creating the chosen 
direction, the held lowers the overall symmetry of the structure what results 
in remarkable phenomena such as, for example, the change in the interband 
optical electroabsorption of a semiconductor QW as compared to the bulk 
samples. This quantum-con£ned Franz-Keldysh effect [8-12] is employed in 
the design of the devices of integrated optics, such as high-speed modulators 
[13], optically bistable switches [14] and wavelength selective detectors [15]. 

The condition 

\k|5 = 0 (1) 


of the zeroing at the confining surface S of the solution T of the wave equation 
is not the only demand implemented in the analysis of the physical and 
chemical systems. Another frequently used type of the boundary condition 
(BC) is the requirement of zeroing at the interface of the normal derivative 
of T: 

(9\k , , 

Here, n is an inward unit normal to the surface. This condition is satis- 
hed, for example, by the transverse electric waves in the electromagnetic 
waveguides [16] or the order parameter of the superconductor bordering the 
vacuum or the insulator [17]. The planar structures with different BCs on 
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the opposite confining walls [18] are also encountered in nature; for example, 
the long-distance dynamics between the Earth and ionosphere of the very low 
frequency electromagnetic oscillations with wavelengths 10 < A < 100 km is 
described, in the first approximation, as a propagation in the fiat waveguide 
with one of its walls (the ground) being a perfect electric conductor with a 
reflection coefficient -1-1 what corresponds to the BC from Eq. (2), while the 
ionosphere behaves as an ideal magnetic conductor with a reflection coeffi¬ 
cient —1 with the edge requirement governed by Eq. (1) [19]. Similarly, the 
sound pressure in the acoustical waveguide formed by the air-water surface 
and the bottom of the lake is zero at the upper free boundary with its van¬ 
ishing derivative at the pond’s rigid floor [20]. The same configuration of 
the BCs is used, e. g., in the description of the continental shelf waves along 
curved coasts [21]. A short review of other structures with mixed combina¬ 
tion of the Dirichlet and Neumann BCs, which are called Zaremba geometries 
[22], is given in Ref. [23]. 

In the present research, exact calculations of the QW with all possible 
permutations of the BCs from Eqs. (1) and (2) in the uniform transverse elec¬ 
tric field are presented and analyzed with the emphasis on the interrelations 
between them. While the case of the Dirichlet, Eq. (1), BCs at the both edges 
is quite well known [2-11], other distributions of the interface requirements 
with a voltage drop between them have received almost no attention [24]. 
Meanwhile, they are indispensable, for example, in the analysis of the pro¬ 
cesses in the metallic waveguide with insulator or ferrite inner filling whose 
refractive index linearly changes between the confining walls. In addition, if 
the electron concentration in the ionosphere is also a linear function of the 
height above the Earth [25], then the corresponding phenomena will be de¬ 
scribed too by the rules presented below. However, our best intention here is 
not only to advice specific applications of the obtained results but, first of all, 
to build up a more general picture, which is also of the large didactic signif¬ 
icance, by establishing interdependencies, drawing parallels and discovering 
similarities and differences between miscellaneous geometries of the BCs and 
their interplay with the field. Following purely Dirichlet case [2-11], we, for 
definiteness, will speak about the quantum particle of the mass m and charge 
—e ( with e being an absolute value of the electron charge) and the wave- 
function that governs its motion through the solution of the corresponding 
steady state Schrodinger equation. A transformation to, e.g.. Maxwell equa¬ 
tions is elementary to do as both these relations belong to the same class 
that in mathematics is called a Helmholtz equation. Excellent discussion 


3 


of the analogies between the Helmholtz and time-independent Schrodinger 
equations [26] shows that drawing such parallels becomes one of the main 
driving ideas in nanophotonics. Upon calculating and analyzing energies 
and functions, a derivation is done of the polarizations of the states, which 
are a quantitative measure of a system response to the held. A fundamental 
relation that directly follows from the Hellmann-Feynman theorem and links 
the polarization with the rate of the energy change with the held, is derived 
and analyzed. Consideration of the asymptotic limits of the weak and strong 
helds allows to arrive at simple analytical expressions. Quantum informa¬ 
tion entropies in position and momentum space are calculated and analyzed 
too; in particular, an andvantage of the entropic uncertainty relation with 
respect to its Heisenberg counterpart is shown for the non-Dirichlet BCs. 
For the zero held, it is proved that the position entropy is BC independent 
and for all states but the lowest Neumann level equals to ln2 — 1 while the 
momentum entropy depends on the level and type of surface requirement. 
The held varies position and momentum entropies in the opposite directions 
in such a way that the corresponding uncertainty relation is always satis- 
hed. Some other results, such as the magnitude of the held at which the 
excited-state polarizations vanish, are discussed too both numerically and 
analytically. Summary of the properties of the eigenvalues and eigenfunc¬ 
tions of the Laplace operator in hat arbitrary bounded Euclidean domains 
with, in particular, Dirichlet, Eq. (1), or Neumann, Eq. (2), BC is given 
in a recent review. Ref. [27]. So, in a sense, the discussion below is a hrst 
extension of the previous research for the simplest structure with diherent 
edge conditions biased by the electric held. A companion paper [28] uses the 
results obtained here for the calculation of the statistical properties of the 
same system. 

The outline below is as follows. In Sec. 2 our physical model is introduced 
and the necessary mathematical formalism for its analysis is developed. Sec. 3 
is devoted to the presentation of the results with their detailed numerical and 
analytical investigation. Sec. 4 wraps up the study by some conclusions. 

2 Model and Formulation 

A charged particle in the ID inhnitely deep QW of the width L is subject to a 
uniform electric held S that is applied perpendicularly to its interfaces. Fig. 1. 
On each of the walls either Dirichlet (D) or Neumann (N) BC is imposed. 
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Accordingly, below we will denote the structure by the two characters where 
the hrst (second) one refers to the BC at the left (right) wall. Of course, ND 
conhguration can be considered as the DN geometry with the electric held 
pointing in the opposite direction. However, to keep the results consistent 
with the uniform BCs that are symmetric with respect to the sign change of 
the held, we discuss positive values of S' only and consider ND and DN cases 
separately. The Hamiltonian H in the Schrodinger equation 

HT(x) = E^iix) (3) 


consists of the held-free part Hq and the electrostatic energy Vs{x) of the 


held 

with 


H = Ho + V4x) 


(P 


Hn = 


2m dx'^ 


+ Vqw{.x) 


(4) 

(5) 


and 


Vs-{x) = -eSx. 


QW potential Vqw{x) reads: 

Vqw{x) 


f 0, |a;| < L/2 

\ oc, |x| > L/2. 


( 6 ) 

(7) 


Our primary goal is to hnd and analyze energies E and wavefunctions \l/(a;) 
for all possible permutations of the BCs and, based on them, to calculate 
other quantities describing the system; in particular, polarization P that is 
dehned as [7] 

P{S) = {ex)^ - {ex)^^Q (8) 

with (...) denoting a quantum mechanical expectation value: 

fL/2 

(x) = / x'^‘^{x)dx. (9) 

J-L/2 

It is convenient from the very beginning to use dimensionless units such 
that the distances are measured in units of the well width L, and energies 
- in units of the ground state energy iPP ji^mL?) of the Dirichlet QW. 
Then, the unit of the electric held iPti}j{2emL^') is naturally determined 
from the voltage drop eSL across the well, and unit of the polarization will 
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Figure 1: Potential profile V{x) of the QW located between —L/2 < x < L/2 
in the uniform electric held S' directed in the negative x direction. Origin of 
the potential axis coincides with the middle of the potential drop across the 
well. At each wall either Dirichlet or Neumann BC is satished. 
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be eL. Hence, the field-free part Hq and electrostatic potential Vg{x) in the 
Hamiltonian H from Eq. (4) are given as 


Hn = 


1 d? 


TT^ dx^ 


+ Vqw{x) 


Vs{x) = -Sx 


with 

Vqw{x) 

BCs are: 

for the Dirichlet QW: 


0, |a:| < 1/2 

oo, |a:| > 1/2. 



for the ND case: 



for the DN case: 



for the Neumann QW: 



( 10 ) 

( 11 ) 

( 12 ) 

(13a) 

(13b) 

(13c) 

(13d) 


with the prime denoting a derivative of the function with respect to its argu¬ 
ment X. Solutions of the Schrodinger equation form a countably infinite set 
with the wavefunctions ^l^{S’]x) for the nth level (n = 0,1,2,...) written 
as the linear combination of the Airy functions Ai(.s) and Bi(z) [29]: 
for the pure Dirichlet BCs: 

H-n - H-n H-n] (i^a) 

for the ND case: 

dBi(.r) Ai(.r) - Ai(.r) Bi(.r)] (^b) 

for the DN case: 

- Ai'(.r) Bi(zr)] (14c) 

for the pure Neumann BCs: 

cr[Bi'(.D Ai(.r) - Ai'(.)r_^) Bi(.r)], (i4d) 
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where 


(15) 

(16) 


.y ^ 4-' KB";-) = -(J)'V 

=±i(>rV)‘'=- (l)'V 

with superscripts I and J taking the values of D and/or N. Expressions for 
the normalization constants C„, which are determined from 

/•1/2 

/ ^l{x)dx = l, (17) 

J-1/2 

will be discussed below. The form of the functions from Eqs. (14) automati¬ 
cally satisfies the BC at the right edge, x = 1/2, for all four configurations. 
Applying the corresponding requirement to the the second conhning surface, 
one gets transcendental equations for finding eigenenergies E^. 
for the pure Dirichlet BCs: 

Ai(e^) - Bi(e^) = 0 (18a) 

for the ND case: 

Ai(e^) Bi'(.„T) - Bi(.r) Ai'(.„T) = 0 (18b) 

for the DN case: 

Ai'(e^) Bi(.„^/) - Bi'(e^) Ai(z„^/) = 0 (18c) 

for the pure Neumann BCs: 

Ai'(.r) Bi'(C^) - Bi'(.r) Ai'(C^) = 0. (18d) 

Eqs. (18) represent corresponding energies {S') as implicit functions E^"^ 
of the electric induction S\ 


E{En, S) — 0. (19) 

Calculation of the coefficients Cn is greatly facilitated by the fact that 
an indefinite integral of the products of the two linear combinations of the 
Airy functions Ai(a: -|- a*) and Bi(x -1- CKj), f = 1, 2, with, in general, different 
constants a\ and is expressed through permutations of the products of the 
same combinations and their derivatives [30]. This, together with the BCs 


from (13) and following from them Eqs. (18) leads to quite simple expressions 
for Cni for example, for the Dirichlet QW it is: 


/^DD _ 

<-n = 



( 20 ) 


where also the value of the Wronskian of the Airy functions, which is equal 
to l/vr [29], has been used. Note that the sign in Eq. (20) is chosen in such 
a way that the functions in the vicinity of the right edge do not become 
negative, \l^„(<#’;x —)■ 1/2) > 0, what is more convenient for our subsequent 
discussion. This completes the calculation of the normalized wavefunctions 
\h„(a;), which correspond to the energies and ushers in the hrst method 
of Ending the polarization that in our dimensionless units reads: 


P(<S) = 


( 21 ) 


It relies on the direct calculation of the integrals in Eq. (21). First, one needs 
to find their values at zero electric field. It is important to underline that 
the well known for this case expressions for the energies 



Pr (0) = (22) 


and functions 



(23a) 


(23b) 


(23c) 


(23d) 


can be derived also with the help of the asymptotic properties of the Airy 


functions [29] as the limit ^ 0 of the corresponding dependencies from 


Eqs. (18) and (14), respectively. Close enough to the right edge, the functions 
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from Eqs. (23) stay positive for any permutation of the BCs and the arbitrary 
n, \l/„(0;a: 1/2) > 0. Elementary calculation yields: 


(x) 


DD _ 
<g’=0 ~ 


{x) 


NN 


0 


(x) 


DN _ 
£’=0 ~ 


(x) 


ND 

£=0 


1 

27r2(n + 1/2)2' 


(24a) 

(24b) 


These expected results are naturally explained by the symmetry of the flat 
QW with respect to its middle for the uniform BCs and the lack of it for the 
Zaremba conhgurations. Applied electric held destroys (for the DD and NN 
structures) this symmetry or changes (for the mixed BCs) the asymmetry 
pushing in this way the polarization from its zero value. Calculation of the 
corresponding integrals in Eq. (21) can again be performed analytically [30] 
producing, after some simple algebra: 


Pn^i^) = ~ 


2Er 


pri'^) = 


pDN^^) ^ 






3 (f 


, lAi^ra+Ai^(^ 
1 


+ 


e 27r2(n +1/2)2 

n/3 + kTBi(e^)/Bi-(.„T)]^ 

1 + z^^ [M{z^F)/M{z^n^ 

1 

(f 27r2(n +1/2)2 

n /3 [B/(.f_^)/Bi(zg/v)]^ 

zr + [Ai'(^r )/Ai(^^/)]^ 

j^NN 

n/3 kTBi'(.r)l'- 


(25a) 


(25b) 


(25c) 


(25d) 


Another way of arriving at the same equations will be given in the next 
section. 

Important characteristics of the physical and chemical systems is their 
quantum information entropy, which can be calculated in the position as 
well as momentum coordinates. For the ID QW under consideration, the 
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( 26 ) 


position and momentum Sk entropies are defined as 

/• 1/2 




'- 1/2 


p{x) In p{x)dx 


Sk = - 


■y{k) \n'y{k)dk, 


(27) 


where p{x) and 7 (/c) are, respectively, position and momentum space densi¬ 
ties: 


p(a;) = |d/(x)p (28) 

^{k) = \m\^ (29) 

with momentum space wave function $(/c) being a Fourier transform of its 
position counterpart: 

m = j e-“'<I<(x)<ix. (30) 

Paramount physical signihcance of these quantities is due to the fact that 
the total information entropy 

St = S^ + Sk (31) 

allows one to express the quantum uncertainty relations in the form 

S'* > 1 -|- In TT. (32) 

Instructive history of the derivation of the d-dimensional analogue of Eq. (32) 
is described in Ref. [31]; namely, the hrst conjectures of these relations from 
the late 1950-ies [32,33] were rigorously conhrmed in 1975 by Bialynicki- 
Birula and Mycielski [34] and Beckner [35]. It is known that the entropic 
uncertainty relation [36, 37] from Eq. (32) is stronger than its Heisenberg 
counterpart that is written in the form 

AxAk > - (33) 

with Ax and Ak being position and momentum standard deviations: 

Ax = \J {pc^) — (x)^ 
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(34) 

(35) 





where the momentum averaging is performed on the basis of the correspond¬ 
ing functions ^{k): 

/ OO 

k^-f{k)dk. (36) 

-CX) 

Shortcoming of Eq. (33) as compared to the entropy relation (32) can be 
shown for our system as well; consider, for example, the held-free lowest 
NN state [36,38] that due to the constant unit magnitude position function 
[x] = 1 exhibits the following momentum density: 

^ 0 sin 0 ixtS = 0. (37) 

Obviously, Heisenberg inequality (33) in this case is meaningless since the mo¬ 
mentum deviation Afc diverges due to the infinite integral for {k'^). However, 
the corresponding entropy is hnite and can be calculated numerically 
producing the result of 2.6834, which, together with the zero position entropy 
= 0, does satisfy inequality (32) with its right-hand side equal to 
2.1447. Note that for Eq. (33) the situation is aggravated even worse if, in¬ 
stead of taking a momentum expectation value on the basis of the functions 
<h(/c) (where A: is a number), one employs the averaging on the basis of the 
position functions 4/(a:) when the momentum is an operator, k = 

/‘V2 

(^)x ^ / '^*{x)k'${x)dx. (38) 

J-l/2 

For the configuration considered above this immediately leads to Ak = 0 
what obviously violates (33). It was shown recently [39] that this discrepancy 
is resolved if one recalls that the Heisenberg inequality (33) was derived for 
the infinite volume where the wave function vanishes at inhnity. For the 
non-Dirichlet edge requirements of the finite ID system a generalized BC- 
dependent uncertainty condition was derived [39] that for the lowest held- 
free Neumann state produces > 0, which is indeed satished as 

an equality. For the Dirichlet conditions, the Heisenberg relation remains 
true since the function \l/(x) vanishes at the surface and, accordingly, at the 
inhnity; moreover, in this case the momentum deviation Ak is the same 
either calculated over the position \[^(a;) or momentum $(/c) functions: 

= + (39) 
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Figure 2: Momentum space densities 7 ^ at zero electric field for (a) the 
Dirichlet, (b) Neumann, and (c) mixed BCs as functions of the momentum 
k where the solid lines are for the ground states (n = 0 ), dashed curves - for 
the first excited levels (u = 1 ), dotted lines - for the second excited states 











Arguments outlined in the previous paragraph dictate a necessity of 
studying quantum information entropies, especially for the structures with 
the non-Dirichlet BCs. Miscellaneous aspects of the properties of the den¬ 
sities p{x) and 7 (/c) and entropies Sx and Sk for different quantum systems 
were addressed in the literature [40-50]. Before discussing the electric field 
influence on the densities and entropies, it does make sense to investigate 
these quantities for the flat QWs. Relevant to the present discussion, let us 
mention that the position and momentum entropies of the Dirichlet struc¬ 
ture were calculated numerically [44,46] and it was shown analytically [43] 
that = ln2 — 1 —0.30685, independent of the quantum number n. 

Applying known integrals from the famous reference book [51], it is easy to 
show that at the zero helds the same is true for all other BCs but the ground 
Neumann state, which, of course, is zero: 

c./j| 0, n = 0,/ = J = iV 

^ 1^=0 In 2 — 1 , all other cases. ^ 

On the other hand, for the momentum entropies at the zero voltage with the 
densities 




ln^{k) 


2(n + l)^ 7 r [1 -f (— 1 )" cos k] 
[/c^ — [n + l)^ 7 r 2 ]^ 

2 /c^ [1 — (—l)"'cosA:] 

TT (/c^ — 


In 


(k) 


2 (—l)"’+i (^ri + nksink + k'^ + [n + tt 


vr 


k^ 


(n+i)V2 


(41a) 

(41b) 


(41c) 


and for any integrals from Eqs. (26), (27) and (30) at (f 7 ^ 0 there are no 
analytical expressions in known to us literature [29,30,51-55]; accordingly, 
their direct numerical quadrature was performed in calculating the results 
presented below. Physically, zero entropy in Eq. (40) means that the uncer¬ 
tainty of determining the position of the particle is the largest one. 

Fig. 2 exhibits zero-field momentum densities juik) from Eqs. (41) and 
(37) for all BCs and several quantum states n. Previous analysis discussed 
these dependencies for the Dirichlet QW only [41,44-46]. As 'Jnik) is an even 
function of its argument, 'jni.—k) = 1 n{k), we plot the parts with the non 
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negative k only. Qualitatively, this behavior is the same for any type of the 
surface requirements; namely, for all excited states, n > 1, the momentum 
density is characterized by the two conspicuous peaks, which are located sym¬ 
metrically with respect to /c = 0 and are flanked by the series of much smaller 
extrema. The distance between the sharp maxima decreases for the smaller 
n, and for the ground state they merge together forming one symmetric 
peak with its magnitude being the largest as compared to the excited lev¬ 
els. Quantitatively, Neumann maximum 7(J^^(0) = l/(27r) = 0.1592 is larger 
than the Dirichlet or Zaremba one, 7(f^(0) = 7(^^(0) = 4/7r^ = 0.1290. As 
the momentum density is smaller than one, 0 < 7 „(/c) < 1, the corresponding 
entropy Sk from Eq. (27) is always positive. More detailed numerical analysis 
reveals that it is a growing function of the quantum number n. In contrast, 
the parts of the position density, which are larger than unity, in their contri¬ 
bution to the integral from Eq. (26) overweigh those with pn{x) < 1 leading in 
this way to the negative entropies from Eq. (40). Returning to the momentum 
entropies, we note that despite the fact that 7 (^^( 0 ) < 7 (^'^( 0 ), the ground- 
state entropy = 2.9000 is greater than its mentioned above Neumann 
counterpart what apparently is explained by the far-reaching spread of the 
corresponding density: the wide sidelobe is clearly seen for the n = 0 curve 
in panel (c) of Fig. 2. In this way, the Dirichlet ground state with its smallest 
momentum entropy = 2.5189 comes closest to the limit imposed by the 
entropic uncertainty relation: = 2.2120. Our calculations indicate 

that the Dirichlet momentum and total entropy for any level n remain the 
smallest ones among all possible BC combinations: . 


3 Results and discussion 

3.1 Energy and polarization 

Energy spectrum as a function of the held S’ is shown in Fig. 3 for all four 
types of the possible combinations of the BCs, and corresponding polariza¬ 
tions for the ground and two lowest excited states are depicted in Fig. 4. It 
is seen that the energies, which at the zero held have the same values for 
the opposite Zaremba geometries, Eq. (22), at the small voltage move in the 
opposite directions while at the large S they coincide with their counter¬ 
parts for the uniform BC distribution that is imposed at the right surface 
of the asymmetric edge requirements. This means that the ND geometry 
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Figure 3: Energy spectrum E as a function of the electric field S’ for the pure 
Dirichlet (solid lines), Neumann (dotted curves), Dirichlet-Neumann (dashed 
lines) and Neumann-Dirichlet (dash-dotted curves) QW. Horizontal straight 
line denotes zero energy. 
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Figure 4: Polarizations Pn for (a) the ground, n = 0, (b) first, n = 1, and (c) 
second, n = 2, excited states as functions of the normalized electric held S . 
Solid (dotted) lines denote pure Dirichfot (Neumann) conhguration while the 
dashed (dash-dotted) curves are for the DN (ND) geometry. Thin horizontal 
lines in parts (b) and (c) denote zero polarizations. Note different P scales 
in each of the panels. 













is the only configuration at which the lowest energy at the small fields in¬ 
creases. Ground-state polarizations for all BCs monotonically increase with 
the field and approach at the strong voltages the asymptotic values that will 
be discussed below while for the excited levels the polarizations take negative 
values for the small and moderate fields, and only after reaching minimum 
they increase with the growing S. Classically, this result, which first was 
predicted for the Dirichlet QW [7], is unexpected since it means that the 
charged particle moves against the applied force. However, in the quantum 
mechanical treatment the electron is described as a wave whose behavior is 
governed by the solution of the corresponding equation inside the domain 
with the appropriate BCs. For the small electric fields, the energy and polar¬ 
ization dependencies are conveniently described by the standard perturbation 
theory [1]; namely, considering in the total Hamiltonian K from Eq. (4) the 
second term as the small disturbance of the field-free Hamiltonian i/o, one 
writes the perturbed energies and wavefunctions as 

oo ' 12 

^ (42, 

m=0 ' ' ' ' 


These expressions are valid at 




•C 1. Here, the prime means that 


^n(O) — Ern{0) 

the term with m = n is dropped out from the series, and are the matrix 
elements of the electric potential, Eq. (11), which are calculated on the basis 
of the undisturbed functions T„(0;a:) from Eqs. (23): 


= 

* mn 


»l/2 


1 - 1/2 


Tm(0; x)Vg{x)'^n{0; x)dx 


/■1/2 

= -S I X'^rn{^]x)'^ri{^]x)dx. 
J-1/2 


(44) 


Note that the condition of the applicability of the perturbation theory re¬ 
duces for the low-lying states basically to ^ 1 that, for simplicity, will 

be used below; however, one has to remember that for the high-lying levels, 
n the approximation produces very good results at the n-dependent 

applied voltage S' that actually can be much larger than unity. Integration 
immediately shows that for the symmetric BCs the energies for the quite 
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small fields are quadratically proportional to S while for the Zaremba ge¬ 
ometries they contain also a linear dependence on the voltage: 


X 


X 




2 , 16(n+l)2 


TT^ 




(m -|- 1)" 


m=0 

ND 


[{n + 1)2 - {m + 1)2]® 


^2 


E^^(<f<l) = (n + ^ ) ±——^ 

V 27 27r2(n + l/2)^ 

1 [(-1)"+™ - 1]^ 

1 ~ m)^(n -|- m -|- 1) 




{n — m){n -|- m -I- 1)® 






2-h 


nO 


TT^ 


°° ('^2 I ™ 2\2 

E(2-^„o) 


m=0 


(45a) 


(45b) 


(45c) 


where Snm “ | g' ^ Kronecker S. Infinite series in Eqs. (45) 

are easily calculated on the basis of the method proposed in Ref. [3] for the 
Dirichlet case. As an example, let us consider its application for the ground 
state, n = 0, of the Zaremba geometry. After some simple algebra, the series 
from (45b) is transformed to 


- «E 


- 128 E 


[(2m)^ - 1]' 
1 


-96E 


l(2m)^ - 1]" 


[(2m.)l - 1] 


5 • 


(46) 


Explicit expression for ~ a^)“^ with the noninteger a is provided 

in Eq. 5.1.25.37 of Ref. [56] with other series being missed there and in other 
sources; however, for our purposes it is useful to point out that all of them 
can be readily calculated by observing that each differentiation with respect 
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to a of the identity [57] 



2a^ 


1 TT 


-cot vra 

2a 


( 47 ) 


produces one additional power of in the denominator of the series 

in its left-hand side. After that, putting a equal to 1/2 allows to hnd the 
expressions above. In the same way, the series ~ are 

found by A; — 1 differentiations of the equality [58] 



(48) 


In addition, the series appearing in the calculation of 

Eq^{S’) can be derived as a limiting case a —)■ 0 of the above mentioned 
summation, or one can use directly the same source [59] to get its value as 
7r®/960. For reference, the expressions for the three lowest levels at -C 1 
are provided below: 



(49a) 


(49b) 


(49c) 


(Eq. (49a) was derived before either with the help of the perturbation calcu¬ 
lations [3, 5] or variational approach [5, 6] or hypervirial-perturbational treat¬ 
ment [4]) 



0.25 + 0.203(f-0.0207(f2 


(50a) 



2.25 + 0.0225<r+0.0104f?2 


(50b) 



6.25 + 0.00811(f + 0.00348(f2 


(50c) 
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(expressions for the DN case are given by Eqs. (50) with, however, the neg¬ 
ative sign of the linear term) 

=-0.0822^2 (51a) 

Ef^(^) = i + + = i + 0.0652^2 

lo \7r"' 8/ 

E^^(ff) = 4 + 4 = 4 + 0.00798^^ (51c) 

64 V 47r^ 3 / 


Note that in this regime the difference Ei — Eq decreases for the ND ge¬ 
ometry and grows with S' for all other BCs what causes different voltage 
dependence of the critical temperatures of the Bose-Einstein systems in the 
corresponding ID box [28]. Invoking regular units, one can say that, for the 
energy separation between the two lowest states, the small electrostatic held 
leads to the decrease of the effective width for all but ND wells. 

For calculating the polarizations P„, one can insert the expressions for the 
wavefunctions from (43) into the general dehnition, Eq. (21), and perform 
the integration with the subsequent evaluation of the series as it just was 
described above. However, there is another applicable at the arbitrary helds 
elegant method of computing them that employs the Hellmann-Feynman 
theorem [60,61] stating: 


dE 

lx 



(52) 


In this equation, A is some continuous parameter of the system. Applying it 
to our conhguration where the role of the variable A is played by the electric 
held one arrives at 


dEn 

dS 


dH 

Is 


(^) . 


(53) 


where the last equality follows straightforwardly from Eqs. (4), (10) and (11). 
In this way, the following expression for the polarization, which is equivalent 
to that given above, Eq. (21), is derived: 


P»(-«’) = -^-Wm- (54) 

This fundamental relation [62] establishes an intimate connection between 
the physical observable quantity (in our case, polarization) and the speed of 
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the energy change (for our geometry, with the electric held). The magnitude 
of the energy itself is a relative value that depends on the choice of the origin 
from which we measure it; for example, for the pure Dirichlet ground state 
it is a continuously increasing function of the held if zero energy coincides 
with the lowest point located at the right edge [7], see Fig. 1. However, the 
diherence between the two energies has an absolute meaning, as Eq. (54) 
vividly manifests. Applying it to Eqs. (45), one gets: 




= 




32(n + l)^ 


TT^ 


oo 


X 


m=0 

(^) = 


(m + 1)" 


[{n + l)^-{m+iyf 
[(_l)n+rn _ 1]2 




1 (^ ~ m)^{n + m + 1) 


m=0 
2 


[(-!)”+"* + 1] 

{n — m){n + m + 1)^ 

2i2-6no) 

TT^ 






m=0 


(n^ — m‘^y 


(55a) 


(55b) 


(55c) 


valid at <C 1. It is important to note that for the uniform BCs the energies 
at the small helds depend quadratically on the voltage while the polariza¬ 
tions are linear functions of S’. This is explained by the symmetry of the QW 
with respect to its middle, which causes that the hrst-order energy perturba¬ 
tion 1/^ vanishes identically. However, the non-diagonal matrix elements of 
V^{x), which enter into the expression for the perturbed function, Eq. (43), 
determine that the latter is deformed and, accordingly, the polarization sur¬ 
vives already in the hrst order of the expansion with respect to the small 
parameter S. For the Zaremba QWs, there is no such symmetry at the zero 
held; accordingly, both energies as well as polarizations are linear functions 
of the voltage. It is also notable from Eq. (55b) that in this regime the po¬ 
larizations of the opposite mixed BCs are equal to each other, as it is also 
seen from Fig. 4 at the small helds where they are hanked by the uniform 
edge requirements. In addition, it is well known that the second-order cor¬ 
rections to the ground-state energy are always negative [1]; as a result, the 
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corresponding polarizations that are defined from Eq. (54) increase. 

The easiest way to explain qualitatively (and to the large degree of pre¬ 
cision, quantitatively) the negative value of the excited-state polarization at 
the small and moderate electric fields lies in the analysis of the first excited 
level for the pure Neumann BCs. In general, the applied voltage mixes all 
zero-field states, as Eq. (43) shows. But for our consideration it suffices to 
take into account the field-induced admixture to Ti produced by the ground 
level since it is the nearest lying state that, accordingly, exerts the largest 
influence: 




NN 


= 2^/2 


(0;a:) + EoW^(0;x) 


f 

2 ^ 

cos TT X- 

- 

L V 

7r2 


< 1 . 


(56) 


This equation manifests that the small electric admixture of the ground state 
subtracts the x-independent term from x). If one recalls that, by our 

convention, this function is nonnegative at the right interface, x = 1/2, 
it means that the applied voltage decreases (increases) the amplitude of 
the right (left) extremum. As a result, the square of the wavefunction 
|^wv (^;^)|2 (bigger) at the right (left) edge of the QW what 

classically means the paradoxical shift of the particle against the applied 
force. However, quantum mechanical treatment considers the electron as 
a wave and the negative polarization is a manifestation of its wave nature 
equivalent to the interference of the modes in the Fabry-Perot resonator [63] 
with gradually changing inside it index of refraction. The decrease of the 
NN function at the right wall for the small fields is clearly seen in the corre¬ 
sponding panel of Fig. 5 that shows \l/i(x) vs. the electric field for all four 
distributions of the BCs. Polarization of the state elementary follows from 
its definition, Eq. (21), and approximate expression for from (56): 


P(^^(^) ^--ff = -0.164^f, < 1. (57) 

Taking into account all higher lying states modifies insignificantly the above 
result, as an application of the Hellmann-Feyman theorem, Eq. (54), to (51b) 
reveals: 

= -\(L + }^S= -0.130^, ^ < 1. (58) 

8 \7r^ 3/ 

Note that the same two-state approximation applied to the ground level 
produces the following expressions for the function 4/^^(x) and polarization 
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pNN. 
^0 ■ 




<S’.j,NN/ 


4 

= 1 H— :tS’ cos7r( X - 

TT^ 


^ < 1 


(59) 


« —g = 0.1642557,?, ,? « 1. (60) 

TT^ 

We retained in the last result more digits than before in order to show its 
excellent coincidence with the one obtained from Eq. (51a), which counts the 
contributions from all states: 




—^ = 0.1644934(f, < 1. 

60 


(61) 


Comparison of Eqs. (60) and (61) proves that for the ground level the domi¬ 
nant admixture comes from the nearest lying excited state while the contri¬ 
butions from the higher levels can be safely neglected. Eq. (59) tells us that 
for the lowest level the application of the held leads to such deformation of 
the particle concentration inside the well that is consistent with the rules of 
the classical mechanics: probability of Ending the electron at the right (left) 
wall increases (decreases) with the electric force. Dependence on the held of 
the ground state wavefunctions '^^^{x) for all possible I and J is shown in 
Fig. 6. 

Returning to the interaction of the right Neumann wall and the held, let 
us note that each of them separately causes the particle to accommodate 
closer to the interface: the surface attracts the electron by creating an ex¬ 
tremum of its wavefunction while the appropriately directed electric force 
simply pushes it closer to the border. However, for any excited state their 
combined inhuence on the particle in some range of the small voltages leads 
to the opposite result, i.e., to its repulsion from the wall. We underline again 
that this is a manifestation of the wave nature of the quantum charge car¬ 
riers. Fig. 5 shows that for all other combinations of the BCs the closest to 
the right wall extremum undergoes the same suppression of its magnitude 
what results in the negative polarizations exhibited by panels (b) and (c) of 
Fig. 4. Physically, their explanation is absolutely identical to the one dis¬ 
cussed above for the purely Neumann QW; namely, it is the held-induced 
admixture of the diherent states that causes their destructive interference. 
Above, we provided the detailed analysis of the NN structure only since it 
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0 - 0.50 - 0.50 50 


Figure 5: Function x) of the first excited level in terms of the distance 

X and electric field S for all possible permutations of the BCs. In each of 
the panels, the corresponding type of the edge requirements is denoted by 
the two characters. To exhibit better the characteristic features of every BC 
configuration, different viewing perspectives are used in the left and right 
parts of the figure. 
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is the simplest from the mathematical point of view, transparent and very 
instructive one. It is also clear why at the small helds the magnitude of the 
polarization for the ground or excited states is the largest for the purely Neu¬ 
mann QW, as Fig. 4 depicts; namely, for the flat structure all extrema have 
the same amplitude equal to 2^/^, see Eqs. (23), and their relative weight 
in the integral from (9) is determined by its distance from the centre of the 
well: the further from the origin, the larger its contribution is. Obviously, 
the held-induced change of the maxima or minima located at the Neumann 
walls will impact the polarization stronger than their Dirichlet counterparts, 
which are repelled from the boundary into the interior of the well. For the 
same reason, Zaremba geometries with one attractive and one repulsive sur¬ 
faces have their polarizations for <C 1 lying in between those two for the 
uniform BCs. 

Further growth of the electric held makes the perturbation theory less and 
less applicable, and Eqs. (18) and (25) should be used for Ending energies 
and polarizations. It is noteworthy to show that the latter dependencies 
are derivable from the former ones with the help of the Hellmann-Feynman 
theorem that holds for the arbitrary S. Namely, for Ending a derivative 
that enters (54) one considers eigenvalue Eqs. (18) as implicit functions of 
the energy E{S) on the electric held with its most general form provided by 
Eq. (19). A theorem of difierentiating implicit functions states [64]: 


dE dE/dg 
Is ~ ~dE/dE' 


(62) 


Applying this rule to each of Eqs. (18), one arrives at the corresponding 
dependence from (25). Thus, there are two equivalent methods of calculating 
the polarization: one implements a direct integration of (54) with the use 
of the analytical properties of the primitives of the Airy functions, and the 
second one relies on the application of the Hellmann-Feynman theorem and 
the rule for Ending a derivative of the implicit functions. 

At the very strong Eelds, it is reasonable to assume that the electron is 
so squeezed to the right surface that it does not ’feel’ the left boundary what 
leads to the problem of the triangular potential with the vertical Dirichlet 
[65] or Neumann right wall. Then, the wavefunction degenerates to 




> 1 . 


(63) 
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Figure 6: The same as in Fig. 5 but for the ground state function 

For the DN geometry, the viewing location is different as compared to other 

BCs. 
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Application of the right BC leads in this limit to the following energies: 


= = an, ^»1 (64a) 

(^) = (^) = - 2 ^^ - ( 7 j «n, ^» 1 , (64b) 

what results in the polarizations 

P ° V ) = 5 + '^»1 ( 65 “) 

TO = i + (65b) 

TO = ^ + (65c) 

n™(.?) = i + ‘^»1- (656) 

Here, a„ and a(j are the nth solutions of equations [29] 

Ai(a„) = 0 (66a) 

Ai'(a(,) = 0, (66b) 


respectively. Their values are provided, for example, in Table 10.13 of 
Ref. [29]. Also, function AiryAiZeros of computer algebra system Maple"'"'^, 
Ref. [66], calculates with any desired precision. Numerically, it is note¬ 
worthy to remark that the absolute value of the difference between the two 
consecutive zeros is smaller for the Airy function in comparison with its 
derivative: 

(^n+1 ^ ( 67 ) 

what leads to the different critical temperatures of the corresponding Bose- 
Einstein condensates [28]. Eqs. (65) show that there are saturation values of 
the polarizations in the asymptotics of the very high helds: 



ND 

Pn^^(oo) 



2 ^ 27r2(n + 1/2)2 ■ 
28 


(68a) 

(68b) 









Due to small power of S’ in (65), this limit is approached quite slowly, as 
Fig. 4 shows. On the contrary, the equal values of the energies for the QWs 
with the same BC at the right surface, Eq. (64), are achieved at relatively 
moderate voltages, which are level dependent: for the higher lying states this 
happens at the larger helds, see Fig. 3. 

Next, let us discuss the electric helds and corresponding energies at which 
the polarizations turn to zero. In general, they should be found from the 
requirement of the vanishing right-hand-side of (25); for example, for the 
purely Dirichlet case, Eq. (25a) yields the following transcendental equation 
for calculating them; 


. (69) 

_L 1 

^ f^ext ~ 

where the superscript at S means that the energy at this electric held reaches 
maximum. As a hrst approximation, the high-held limit from Eq. (65a) 
provides the value of this voltage as 


^ext 



(70) 


while the corresponding maximum energy at this point is 


rpext 


1 

3 



(71) 


where, for convenience, we have substituted the superscript denoting a type 
of the BC by the one denoting that this is an extremum value. One gets from 
Eq. (70): S^^^ 16.407 and S^^^ 40.408 while the exact ones are: 13.370 

and 38.700. The corresponding approximate energies read: 4.102 and 

El** 10.102 with the exact ones being 4.212 and 10.146. So, the precision 
of the approximation increases with the level number and is better for the 
energy than for the held. It directly follows from Eqs. (70) and (71) that 


4 


rpext 

■^n 

^ext 


1 = 0 , 


(72) 


what automatically zeroes the right-hand side of Eq. (69), as it should be, 
since it contains a contribution from the left wall that is neglected in our 
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high-field approach developed above. Assuming that the expression from 
Eq. (72) is a small nonzero number, 


- 1 = ^, |5| « 1, (73) 

and applying the corresponding Taylor expansion of the left-hand side of 
Eq. (69) with respect to the small parameter 6, one can improve the approx¬ 
imations from Eqs. (70) and (71); however, the obtained results are not very 
transparent and, accordingly, not written here. 

As a next example, it is instructive to introduce zero-energy electric field 
i.e., the field at which the energy of the nth level turns to zero. Its phys¬ 
ical meaning is obvious from its definition and the transcendental equation 
for its determination follows from Eqs. (18) by putting = 0. Assuming 
that this happens at quite large fields, one can invoke asymptotic properties 
of the Airy functions [29] to get 


/a(0) 

®n+l 


DD 


NN 



(74a) 

(74b) 


with the obvious requirement = 0. Even though these dependencies 

NN 

are derived for the high voltages, they produce reasonably good results even 
for the small n; for example, the lowest DN level reaches zero at = 1.118 
while Eq. (74b) calculates it as 1.125. 


3.2 Entropies 

Ground state momentum density is shown in Fig. 7 as a function of the 
momentum k and voltage S'. The general tendency of the electric field in¬ 
fluence is the decrease of the maximal values of 7 „ and smoothing out of the 
oscillations along the k axis. Note that for the Neumann-Dirichlet well the 
largest value 7o(0) as a function of the applied voltage has a local maximum 
of 0.1448 at ^ ~ 3 while for all other geometries it is a monotonically de¬ 
creasing function of the field. Similar dependencies are characteristic for the 
higher lying states too. 
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Figure 7: Ground state momentum density 70 k) in terms of the momen¬ 
tum k and electric field S’. 
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Figure 8: Position Sx (dotted lines), momentum Sk (dashed lines) and total 
St (solid curves) entropies as a function of the electric held S’ for all permu¬ 
tations of the Dirichlet and Neumann BCs. Two characters in the left lower 
corner of each panel denote a type of the edge requirements while the digits 
near the curves depict the corresponding quantum number n. Note vertical 
line breaks from S = 0.01 to S' = 2.19. 
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Peculiar features of the position and momentum densities exemplified by 
Figs. 5-7 determine the properties of the corresponding entropies and 
Sk- Fig. 8 shows them together with the sum St as functions of the electric 
field for all four possible BCs and several low-lying levels. Characteristic 
feature is a nonmonotonic dependence of the position Sx and momentum Sk 
entropies on the small and moderate voltages for the structures with the left 
Neumann wall (except the ground NN level). For example, the maximum of 
the lowest ND state position entropy is explained by the sharp decrease of the 
corresponding function \I'o(2^) at the left surface for < 5 (see the left lower 
panel of Fig. 6) and accompanying decline of its negative contribution to Sx 
while the corresponding density at the right edge stays almost unchanged. 
Simultaneously, mentioned above increase of the momentum density results 
in the contraction of Sk and, as a consequence of these two facts, the lowering 
of the total entropy that, however, stays above the fundamental limit from 
the right-hand side of Eq. (32), as expected. For the higher voltages, the 
accumulating near the right interface function \ho pushes the position space 
entropy Sx downward forming the extremum seen in the figure. On the 
contrary, the Dirichlet left interface subdues the entropy swaying, as the 
corresponding panels of Fig. 8 demonstrate. At the large S’, the entropies 
are determined by the right-wall BC, as discussed above; accordingly, they are 
almost equal for the DN and pure Neumann QWs and their growing with 
the held magnitudes are larger than their (almost identical) counterparts 
for the ND and pure Dirichlet structures. Similar to the held-free case, the 
Dirichlet ground state has the lowest total entropy among all other levels and 
all possible BC combinations. Note its quite small variation with the voltage: 
from its zero-held value of 2.2120 it saturates to approximately 2.2552 at the 
large S. Growing with the held magnitude of the negative position entropy 
Sx means smaller uncertainty in determining particle location while the larger 
values of the positive momentum entropy Sk indicate that the uncertainty in 
calculating momentum k increases. 

4 Concluding remarks 

Exact solutions of the Helmholtz equation inside the ID interval with dif¬ 
ferent BCs at its ends and subjected to the linear potential inside revealed 
similarities as well as diherences of the miscellaneous permutations of the 
edge requirements. Similarities include: a continuous growth with the elec- 
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trie field S’ of the ground-state polarization for any BC geometry and its 
slow approach to the saturation value that, however, is different for the uni¬ 
form and Zaremba distributions of the surface requirements; a decrease of 
the excited-level polarization for the small voltages and reverse of its change 
for the stronger electric forces. The negative polarization at the small S is 
conveniently described by the perturbation theory; in particular, it is shown 
that its reason lies in the held-induced mixing of the unperturbed levels. 
Mathematical analysis of the simple two-level approximation proved that 
only an admixture of the neighbouring states plays a dominant role in the 
change of the polarization. In addition, the perturbation theory conhrms 
mathematically a physical conclusion about different magnitudes of the po¬ 
larization for the different BCs. In the opposite limit of the high voltages, 
physically the quantum particle inside the box does not ’see’ one of the sur¬ 
faces what mathematically means that solutions of the Schrodinger equation 
in this case depend only on the BC at the opposite wall to which it is pushed 
by the electric held. Application of the Hellmann-Feynman theorem, which 
long ago proved to be a powerful tool of studying quantum systems [60,61], 
in our case revealed a fundamental relation between the speed of the en¬ 
ergy change with the held, on the one hand, and the polarization, on the 
other. Applied voltage also modihes position and momentum densities and 
entropies in such a way that the entropic uncertainty relation, Eq. (32), is 
always satished. 

Dirichlet, Eq. (1), and Neumann, Eq. (2), conditions are the limiting cases 
of the so called Robin BC [67] 

nVT 

where the variation of the extrapolation length A continuously changes the 
edge requirements from the Dirichlet (A = 0) to the Neumann (A = oo) 
conhguration. At the zero voltage, S = 0, the vanishingly small negative 
Robin length, A —)■ —0, heavily pushes two levels to the QW interfaces with 
their energies unrestrictedly decreasing: Ei ^2 — -1 —l/(7rA)^ [39,68]. Their 

almost full degeneracy in this regime might lead to the nontrivial dependence 
of the spectrum on the held while the strong localization at the QW walls 
might produce unusual polarization and entropy behavior. As a simplifying 
example, for the single wall with the negative Robin length A [69-72] it is 
elementary to prove that at = 0 it possesses one bound state with the 




(75) 
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energy (in regular units) 


2mA2 


(76) 


E = 


and the entropies 


Sx = 1 — In 2 + In |A| 

Sk = 21n2 + luTT — In |A|, 


(77a) 

(77b) 


where it is assumed that ff/{2m) = 1. Eqs. (76) and (77) are exact for 
any negative A. It is seen from them that even though the position and 
momentum entropies diverge at the limiting values of A ^ — cxo and A —)■ —0, 
the total entropy 

S'* = 1 + luTT + ln2 (78) 

stays hnite, does not depend on the extrapolation length (due to the exact 
cancelling of the divergent terms in Sx and Sk) and does satisfy inequality 
(32). 

At the weak helds the polarization for any type of the BCs has opposite 
signs for the ground and any excited states. The natural question arises: Is 
there “warm” enough temperature or/and sufficient number of the charged 
fermions that are able to push the total statistically averaged polarization to 
the negative values? The following paper [28], among other issues, answers 
this question. 
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